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Abstract 

We consider a financial market with a stock exposed to a counterparty risk inducing 
a drop in the price, and which can still be traded after this default time. We use a 
default-density modeling approach, and address in this incomplete market context the 
expected utility maximization from terminal wealth. We show how this problem can 
be suitably decomposed in two optimization problems in complete market framework: 
an after-default utility maximization and a global before-dcfault optimization prob- 
lem involving the former one. These two optimization problems are solved explicitly 
respectively by duality and dynamic programming approaches, and provide a fine un- 
derstanding of the optimal strategy. We give some numerical results illustrating the 
impact of counterparty risk and the loss given default on optimal trading strategies, in 
particular with respect to the Merton portfolio selection problem. 



Key words: Counterparty risk, density of default time, optimal investment, duality, dy- 
namic programming, backward stochastic differential equation. 
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1 Introduction 



In a financial market, the default of a firm has usually important influences on the other 
ones. This has been shown clearly by several recent default events. The impact of a coun- 
terparty default may arise in various contexts. In terms of credit spreads, one observes 
in general a positive "jump" of the default intensity, called the contagious jump and in- 
vestigated firstly by Jarrow and Yu [5]. In terms of asset (or stock) values for a firm, the 
default of a counterparty will in general induce a drop of its value process. In this paper, 
we analyze the impact of this risk on the optimal investment problem. More precisely, we 
consider an agent, who invests in a risky asset exposed to a counterparty risk, and we are 
interested in the optimal trading strategy and the value function when taking into account 
the possibility of default of a counterparty, together with the instantaneous loss of the asset 
at the default time. 

The global market information containing default is modeled by the progressive enlarge- 
ment of a background filtration, denoted by F, representing the default-free information. 
The default time r is in general a totally inaccessible stopping time with respect to the 
enlarged filtration, but is not an F-stopping time. We shall work with a density hypothesis 
of the conditional law of default given F. This hypothesis has been introduced by Jacod [4] 
in the initial enlargement of nitrations, and has been adopted recently by El Karoui et al. 
[3] in the progressive enlargement setting for the credit risk analysis. The density approach 
is particularly suitable to study what goes on after the default, i.e., on {r < t}. For the 
bef ore-default analysis on {r > t}, there exists an explicit relationship between the density 
approach and the widely used intensity approach. 

The market model considered here is incomplete due to the jump induced by the default 
time. The general optimal investment problem in an incomplete market has been studied 
by Kramkov and Schachermayer [7] by duality methods. Recently, Lim and Quenez [8] 
addressed, by using dynamic programming, the utility maximization in a market with 
default. The key idea of our paper is to derive, by relying on the conditional density 
approach of default, a natural separation of the initial optimization problem into an after- 
default one and a global before-default one. Both problems are reduced to a complete 
market setting, and the solution of the latter one depends on the solution of the former 
one. These two optimization problems are solved by duality and dynamic programming 
approaches, and the main advantage is to give a better insight, and more explicit results 
than the incomplete market framework. The interesting feature of our decomposition is to 
provide a nice interpretation of optimal strategy switching at the default time r. Moreover, 
the explicit solution (for the CRRA utility function) makes clear the roles played by the 
default time r and the loss given default in the investment strategy, as shown by some 
numerical examples. 

The outline of this article is organized as follows. In Section 2, we present the model 
and the investment problem, and introduce the default density hypothesis. We then explain 
in Section 3 how to decompose the optimal investment problem into the before-default and 
after-default ones. We solve these two optimization problems in Section 4, by using the 
duality approach for the after-default one and the dynamic programming approach for the 
global before-default one. We examine more in detail the popular case of CRRA utility 
function and finally, numerical results illustrate the impact of counterparty risk on optimal 
trading strategies, in particular with respect to the classical Merton portfolio selection 
problem. 
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2 The conditional density model for counterparty risk 



We consider a financial market model with a riskless bond assumed for simplicity equal to 
one, and a stock subject to a counterparty risk: the dynamics of the risky asset is affected 
by another firm, the counterparty, which may default, inducing consequently a drop in 
the asset price. However, this stock still exists and can be traded after the default of the 
counterparty. 

Let us fix a probability space (Q, Q, P) equipped with a brownian motion W = (Wt)t£\o,T\ 
over a finite horizon T < oo, and denote by F = (J-t)te[o,T\ the natural filtration of W. We 
are given a nonnegative and finite random variable r, representing the default time, on 
(fi, £7,1P). Before the default time r, the filtration F represents the information accessible 
to the investors. When the default occurs, the investors add this new information r to the 
reference filtration F. We then introduce Dt = l T <t, < t < T, D = {T^t)te[o,T] the filtration 
generated by this jump process, and G = (Gt)te[o,T\ the enlarged progressive filtration FVB, 
representing the structure of information available for the investors over [0, T\. 

The stock price process is governed by the dynamics: 

dS t = S t -(titdt + a t dWt-7tdDt), < t < T, (2.1) 

where fi, a and 7 are G-predictable processes. At this stage, without any further condition 
on the default time r, we do not know yet that W is a G-semimartingale (see Remark 12. ip . 
and the meaning of the sde (12. ip is the following. Recall (cf. [9]) that any G-predictable 
process <p can be written in the form: tpt = ipJlt< T + (pf(r)lt> T , < t < T, where ip v is 
F-adapted, and <pf(9) is measurable w.r.t. Tt <S> £>(M+), for all t € [0,T]. The dynamics 
(|2.ip is then written as: 

dS t = S t (f4dt + afdW t ), 0<t<r, (2.2) 
S T = 5 r -(l- 7 f), (2.3) 
dS t = S t (rf{T)dt + a d t {T)dW t ) : r<t<T, (2.4) 

where cr F , 7 F are F-adapted processes, and (uj,6) -» (4(0), af{9) are T t <g> B([0, t))- 
measurables functions for all t E [0,7"]. The nonnegative process 7 represents the (pro- 
portional) loss on the stock price induced by the default of the counterparty, and we may 
assume that 7 is a stopped process at r, i.e. 7 t = jtAr- By misuse of notation, we shall 
thus identify 7 in (|2.ip with the F-adapted process 7 F in (|2.3|) . When the counterparty 
defaults, the drift and diffusion coefficients (/i, cr) of the stock price switch from (/i F ,cr F ) 
to (fi d (r) , a d (r)) , and the after-default coefficients may depend on the default time r. We 
assume that at > 0, < t < T, the following integrability conditions are satisfied: 
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<Jt| < 00, a.s. (2-5) 

7o °t Jo 
and 

< 7t < 1, < t < T, a.s. (2.6) 

which ensure that the dynamics of the asset price process is well-defined, and the stock 

price remains (strictly) positive over [0,T] (once the initial stock price So > 0), and locally 
bounded. 
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Consider now an investor who can trade continuously in this financial market by holding 
a positive wealth at any time. This is mathematically quantified by a G-predictable process 
7r = (TTt)te[o,T]i called trading strategy and representing the proportion of wealth invested 
in the stock, and the associated wealth process X with dynamics: 

J Q 

dX t = nX t --±, 0<t<T. (2.7) 

By writing the G-predictable process it in the form: 7Tj = 7rflj< r + 7r^(r)lf >T , < t < T, 
where vr F is F-adapted and 7if (0) is Tt®B([0, immeasurable, and in view of (pT2]> - ([2T3]) - <|2~i|> . 
the wealth process evolves as 

dX t = X t 7rf(i4dt + afdWt), < t < r, (2.8) 
X T = X t -(1-4j t ) (2.9) 
dX t = Xf4(T)(t4(T)dt + af(T)dW t ), r<t<T. (2.10) 

We say that a trading strategy it is admissible, and we denote ir E A, if 

T 

{'KtO't^dt < oo, and tt t j t < 1 a.s. 

This ensures that the dynamics of the wealth process is well-defined with a positive wealth 
at any time (once starting from a positive initial capital Xq > 0). 

In the sequel, we shall make the standing assumption, called density hypothesis, on the 
default time of the counterparty. For any t € [0, T], the conditional distribution of r given 
Tt admits a density with respect to the Lebesgue measure, i.e. there exists a family of 
Tt ® £>(]R + )-measurable positive functions (oj, 0) — > Ut(0) such that: 

(DH) F[r£de\Ft] = a t (6)d6, t e [0,T]. 

We note that for any 9 > 0, the process {att(0),0 <t< T} is a (P, F)-martingale. 

Remark 2.1 Such a hypothesis is usual in the theory of initial enlargement of filtration, 
and was introduced by Jacod The (DH) Hypothesis was recently adopted by El Karoui 
et al. [3] in the progressive enlargement of filtration for credit risk modeling. Notice that 
in the particular case where the family of densities satisfies ot(^) = <xt{p) f° r all < t < T, 
we have P[r > t\J~t] = P[r > t|^r]- This corresponds to the so-called immersion hypothesis 
(or the H- hypothesis), which is a familiar condition in credit risk analysis, and means 
equivalently that any square-integrable F-martingale is a square-integrable G-martingale. 
The H-hypothesis appears natural for the analysis on before-default events when t < r, 
but is actually restrictive when it concerns after-default events on {t > r}, see [3] for a 
more detailed discussion. By considering here the whole family {ott(0),t € [0, T],9 € 
we obtain additional information for the analysis of after-default events, which is crucial 
for our purpose. 

Let us also mention that the classical intensity of default can be expressed in an explicit 
way by means of the density. Indeed, the (P, G)-predictable compensator of D t = l T <t is 
given by f^ Ar ag(6)/God9, where Gt = P[r > t\Ft] is the conditional survival probability. 
In other words, the process M t = D t — Jq At aQ(9)/Ggd9 is a (P, G)-martingale. Thus, by 
observing from the martingale property of {at(9),0 < t < T} that Gt = J t at(9)d9 = 
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J" t °° E[ag(9)\Ft]dO, we recover completely the intensity process Xf = lt< T a t (t)/Gt from the 
knowledge of the process {at(t),t > 0}. However, given the intensity A G , we can only 
obtain some part of the density family, namely at{0) for 6 > t. 

Under (DH) Hypothesis, a (P, F)-brownian motion W is a G-semimartingale and admits 
an explicit decomposition in terms of the density a given by (see [9] , [6] , [3] ) : 

Wi _ #o + r« + M = :W? + A„ 0<t<T, 

Jo Jr OL s (t) 

where W G is a (P, G)-brownian motion, and A is a finite variation G-adapted process. 
Moreover, by the Ito martingale representation theorem for brownian filtration F, At is 
written in the form A t = J * a s ds for some G-adapted process a = (at)te[o,T]- Let us then 
define the G-adapted process 

0t = ^ + g ^-^ G , 0<t<T, 

and consider the Doleans-Dade exponential local martingale: Zf = £{— f (3dW G )t, < t < 
T. By assuming that Z G is a (P, G)-martingale (which is satisfied e.g. under the Novikov 
criterion: E[exp(/ ±|/3 f | 2 cft)] < oo), this defines a probability measure Q equivalent to P 
on Qt) with Radon-Nikodym density: 

§ - ^ - -p(-jf -ijflAP*). 

under which, by Girsanov's theorem (see [I] Ch.5.2), W = W G + J /3dt is a (Q,G)- 
Brownian motion, M is a (Q, G)-martingale, so that the dynamics of 5 follows a (Q, G)-local 
martingale: 

dS t = S t -(a t dWf -jtdMt). 
We thus have the "no-arbitrage" condition 

M(G) := {Q ~ P on (n, Qt) : S is a (Q, G) - local martingale} / 0. (2.11) 

3 Decomposition of the utility maximization problem 

We are given an utility function U defined on (0, oo), strictly increasing, strictly concave 
and C 1 on (0, oo), and satisfying the Inada conditions U'(0 + ) = oo, [/'(oo) = 0. The 
performance of an admissible trading strategy ir € A associated to a wealth process X 
solution to (12. 7p and starting at time from Xq > 0, is measured over the finite horizon T 
by: 

Jo(tt) = E[U(X T )}, 
and the optimal investment problem is formulated as: 

Vq = supJ (7r). (3.1) 

■k&A 
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Problem (|3.1|) is a maximization problem of expected utility from terminal wealth in an 
incomplete market due to the jump of the risky asset. This optimization problem can be 
studied by convex duality methods. Actually, under the condition that 



Vo < oo 



(3.2) 



which is satisfied under (|2.1ip once 



^ U(y—-) < oo, for some y > 0, 
L air ' J 

where U (y) = sup x> o [U (x) — xy] , and under the so-called condition of reasonable asymptotic 
elasticity: 



AE(U) := limsup 



xU'{x) 
U{x) 



< 1 



we know from the general results of Kramkov and Schachermayer [7] that there exists a 
solution to (|3.ip . We also have a dual characterization of the solution, but this does not 
lead to explicit results due to the incompleteness of the market, i.e. the infinite cardinality 
of Ai (G) . One can also deal with problem (|3.1|) by dynamic programming methods as done 
recently in Lim and Quenez [8] under (H) hypothesis, but again, except for the logarithmic 
utility function, this does not yield explicit characterization of the optimal strategy. We 
provide here an alternative approach by making use of the specific feature of the jump 
of the stock induced by the default time under the density hypothesis. The main idea is 
to separate the problem in two portfolio optimization problems in complete markets: the 
after-default and before-default maximization problems. This gives a better understanding 
of the optimal strategy and allows us to derive explicit results in some particular cases of 
interest. 

The derivation starts as follows. First notice that any ir G A, thus in the form: 7Tj = 
vrf lt< T + 7T t( T )^t>T, < t < T, can be identified with a pair (tt v , 7r d ) G Aw X Ad where ^4f 
is the set of admissible trading strategies in absence of defaults, i.e. the set of F-adapted 
processes 7r F s.t. 



and Ad is the set of admissible trading strategies after default at time r = 9, i.e. the set of 
Borel family of F-adapted processes {7if(0),0 < * < T} parametrized by 6 G [0,T] s.t. 



Hence, for any ir = (vr F , ir d ) G A, we observe by ([278]) - (f^9]) - (f27T0]) that the terminal wealth 
is written as: 




and T^ele < 1> < 9 < T, a.s. 



(3.3) 




< oo 



a.s. 



where X 



is the wealth process in absence of default, governed by: 



dX: 



■F 



Xfnf (fjfidt + afdWt) , < t < T. 



(3.4) 
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starting from X$ = X , and {Xf(9),9 < t < T} is the wealth process after default occuring 
at r = 9 G [0, T], governed by: 



dX d {9) = X?(9)4(9)(rf(9)dt + at(9)dW t ), 9<t<T 
X d {9) = X ¥ (l-n v j e ). 



(3.5) 
(3.6) 



Therefore, under the density hypothesis (DH), and by the law of iterated conditional 
expectations, the performance measure may be written as: 



Jo(tt) = E[E[U(X T )\F T ]] = E[U{X ¥ )¥[t>T\T t ]+E[U(X^{t))I t < t \T t ]] 



E 



U{Xl)G T + I U(X$(9))a T (9)d9 



(3.7) 



where G T = P[r > T\T T ) = f™ a T (9)d9. 

Let us introduce the value-function process of the "after-default" optimization problem: 



V e d (x) 

4(x,7T d (9)) 



esssup Ji(x,ir d (9)), (9, x) 6 [0, T] X (0, oo) 



(3.« 



E[U(X«> x (9))a T (9)\F e ], 



where Ad{9) is the set of (jF t ) 0<t < T -adapted processes {irf(9),9 < t < T} satisfying 
jj \-K d (9)af{9)\ 2 dt < oo a.s., and {xf ,x {9), 9 < t < T} is the solution to <^ controlled 
by TT d {9) G Ad{9), starting from x at time 9. Thus, V d is the value-function process of an 
optimal investment problem in a market model after default. Notice that the coefficients 
(fi d ,a d ) of the model depend on the initial time 9 when the maximization is performed, 
and the utility function in the criterion is weighted by ax{9). We shall see in the next 
section how to deal with these peculiarities for solving f)3.8[) and proving the existence and 
characterization of an optimal strategy. 

The main result of this section is to show that the original problem (|3. 1 j) can be split 
into the above after-default optimization problem, and a global optimization problem in a 
before-default market. 

Theorem 3.1 Assume that Vg(x) < oo a.s. for all (9,x) G [0, T] X (0, oo). Then, we have: 



In 



sup E 



T 



U{X ¥ T )G T + V e d {X ¥ e {l-4 le ))d9 



(3.9) 



Proof. Given ir = (-/r F , ir d ) G A, we have the relation (|3.7[) for Jo(vr) under (DH). Further- 
more, by Fubini's theorem, the law of iterated conditional expectations, we then obtain: 



Jn(vr) 



= E 
= E 
< E 



U(X$)Gt+ / E[U(X$(9))a T (9)\F e ]d9 



T 



U(X ¥ T )G T + / J d (X d (9)^ d (9))d9 



(3.10) 



U{Xl)G T + 



< sup E 



V d (X d (9))d9 

'0 



U{Xt)G t + I V e d (X ¥ e {l-4le))d9 
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by definitions of J d , V d and X d {9). This proves the inequality: Vo < Vq. 

To prove the converse inequality, fix an arbitrary 7r F G Aw- By definition of V d , for any 
u) G U, 9 G [0,T], and e > 0, there exists Tr d ' £,UJ (9) G «4.d(^)> which is an e-optimal control 
for V^ rf at (uj , Xg (ui , 9)) . By a measurable selection result (see e.g. |10j). one can find 7r d ' e 
G „4 d s.t. vr d ' e (w, 0) = 7r d,B ' w (w, 9), dF (g d0 a.e., and so 

< J${X$(9),it d ' £ (9)), dF®d9a.e. 

By denoting 7r £ = (7r F , 7r d,£ ) G A, and using again (|3.10p . we then get: 

r r T 
Vo > J (vr £ ) = E U(X ¥ T )G T + / 4(X d (9), ^ £ (9))d9 
1 Jo 

~U(X V T )G T + [ V d (X d (9))d9 
1 Jo 



> E 



e. 



From the arbitrariness of tt f in Aw and e > 0, we obtain the required inequality and so the 
result. □ 

Remark 3.1 The relation (|3.9p can be viewed as a dynamic programming type relation. 
Indeed, as in dynamic programming principle (DPP), we look for a relation on the value 
function by varying the initial states. However, instead of taking two consecutive dates 
as in the usual DPP, the original feature here is to derive the equation by considering 
the value function between the initial time and the default time conditionnally on the 
terminal information, leading to the introduction of an "after-default" and a global before- 
default optimization problem, the latter involving the former. Each of these optimization 
problems are performed in complete market models driven by the brownian motion and 
with coefficients adapted with respect to the brownian filtration. The main advantage of 
this approach is then to reduce the problem to the resolution of two optimization problems 
in complete markets, which are simpler to deal with, and give more explicit results than the 
incomplete market framework studied by the "classical" dynamic programming approach 
or the convex duality method. 

Furthermore, a careful look at the arguments for deriving the relation (|3,9p shows that 
in the decomposition of the optimal trading strategy for the original problem (|3.ip which 
is known to exist a priori under (|2.1ip : 

TT t = 7rflt< T +7T d {T)l t>T , 0<t<T, 

7r F is an optimal control to (|3.9|) . and Tt d (9) is an optimal control to Vg(Xg (9)) with Xg(9) 
= Xq(1 — 7r F 7£)), and X ¥ is the wealth process governed by 7r F . In other words, the optimal 
trading strategy is to follow the trading strategy 7r F before default time r, and then to 
change to the after-default trading strategy 7r d (r), which depends on the time where default 
occurs. In the next section, we focus on the resolution of these two optimization problems. 



4 Solution to the optimal investment problem 

In this section, we focus on the resolution of the two optimization problems arising from the 
decomposition of the initial utility maximization problem. We first study the after-default 
optimal invesment problem, and then the global before-default optimization problem. 
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4.1 The after-default utility maximization problem 



Problem (|3.8p is an optimal investment problem in a complete market model after default. 
A specific feature of this model is the dependence of the coefficients (fi d , a d ) on the initial 
time 9 when the maximization is performed. This makes the optimization problem time- 
inconsistent, and the classical dynamic programming method can not be applied. Another 
peculiarity in the criterion is the presence of the density term ot(9) weighting the utility 
function U. 

We adapt the convex duality method for solving (|3.8|) . We have to extend this martingale 
method (in complete market) in a dynamic framework, since we want to compute the value- 
function process at any time 9 £ [0, T]. Let us denote by: 



Zt(9) 



exp 



<{0) 



/4W 



<{0) 



du 



<t<T, 



the (local) martingale density in the market model (|2.4j) after default. We assume that for 
all 9 £ [0, T], there exists some yg ^-measurable strictly positive random variable s.t. 



E 



To 



< oo. 



(4.1) 



This assumption is similar to the one imposed in the classical (static) convex duality method 
for ensuring that the dual problem is well-defined and finite. 



Theorem 4.1 Assume that (|4.ip and AE(U) < 1 hold true. Then, the value- function 
process to problem (13.80 is finite a.s. and given by 

V d (x) = E[u(l{y e (x)^^jy T (e)\^g\, (9, x) € [0, T] x (0, oo), 



and the corresponding optimal wealth process is equal to: 



X?> x (9) 



E 



Zt{9) 



a T {9) 



Ft 



<t<T, 



(4.2) 



where I = (U') 1 is the inverse of U' , and yg(x) is the strictly positive J-q ® £>((0, oo))- 
measurable random variable solution to Xq' x {9) = x. 

Proof. First observe, similarly as in Theorem 2.2 in [7], that under AE(U) < 1, the 
validity of (|4.ip for some or for all yg J-e-measurable strictly positive random variable, is 
equivalent. By definition of Z(9) and Ito's formula, the process {Zt(9)X t ' x (9), 9 <t< T} 
is a nonnegative (P, (jF t )g< t <T)-local martingale, hence a supermartingale, for any ir d (0) 
€ Ad(9), and so E[X^ X (9)Z T {9)\T e ] < xf x {9)Z e {9) = x. Denote Y T {9) = Z T (9)/a T {9). 
Then, by definition of U, we have for all yg J^-measurable strictly positive random variable, 
and ir d (9) G A d {9): 



E[U(X d ' x (9))a T (9)\Fg] 



< 



< 



E[U(y e Y T (9))a T (9)\rg]+E[X d ' x (9)ygY T (9)a T (9)\T e ] 
E[U(y e Y T (9))a T (9)\Fe] + y e E[X d ' x (9)Z T (9)\^ e } 
E[U(y e Y T (9))a T (9)\Te} + xy e , (4.3) 
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which proves in particular that Vg(x) is finite a.s. Now, we recall that under the Inada 
conditions, the supremum in the definition of U(y) is attained at I, i.e. U(y) = U(I(y)) — 
yl(y). Prom (|4.3p . this implies 

E[U(X d ' x (9))a T (9)\T } < E[U{I{y e Y T {9)))a T (9)\T e ] 

- ye(m[Z T (0)I{yoY T (9))\F e }-x). (4.4) 



Now, under the Inada conditions, (|4.ip and AE(U) < 1, for any u) G 0, 9 G [0,T], the 
function y G (0, oo) — ► f$(uo,y) = E[ZT(9)I(yYr(6))\J r e] is a strictly decreasing one-to-one 
continuous function from (0, oo) into (0, oo). Hence, there exists a unique yg(to, x) > s.t. 
fo(u,y(x)) = x. Moreover, since f$(y) is Tq ® 13(0, oo)-measurable, this value ye(x) can be 
chosen, by a measurable selection argument, as J-g (8> Z3(0, oo)-measurable. With this choice 
of ye = ye{x), and by setting X^ x (9) = I(ye(x)Yr(9)), the inequality (|4.4jl yields: 

E[[/(X^(0))a T (fl)|^] < E[[/(Xj*(0))a T (0)|j^], Vtt>) G A^S). (4.5) 

Consider now the process X d ' x (9) defined in (jMJ) leading to X% x (9) at time T. By defi- 
nition, the process {M t {9) = Z t (9)X d ' x (9),9 < t < T} is a strictly positive (P, (J" t ) e < t <T)- 
martingale. From the martingale representation theorem for brownian motion filtration, 
there exists an (^t)e<t<T-adapted process (4>t)e<t<T satisfying L \(f>t\ 2 dt < oo a.s., and 
such that 

M t {9) = M e {9) + [ (j> u M u {9)dW u , 9<t<T. 
Je 

Thus, by setting fc d (9) = (4>+ ^j)/cr d (9), we see that n d (9) G A d (9), and by Ito's formula, 

X d > x (9) = M(9)/Z(9) satisfies the wealth equation f|3.5|) controlled by 7r d (9). Moreover, by 
construction of ye(x), we have: 



xf x (9) = E\Z T {9)l{y e {x) 



Z T {9) \ 
l a T {9)J 



To 



Recalling (|4.5|) . this proves that TT d {9) is an optimal solution to (|3.8|) . with corresponding 
optimal wealth process X d,x (9). □ 

Remark 4.1 Under the (H) hypothesis, ax (9) = ctg(9) is ^/-measurable. In this case, 
the optimal wealth process to (|3.8p is given by: 



X d t' X {0) = E[|^/(w(x)Zr(fl)) 



< t < T, 



where is the strictly positive Tq <g) £>((0, oo)) -measurable random variable satisfying 

Xq' x {9) = x. Hence, the optimal strategy after-default does not depend on the density of 
the default time. 

We illustrate the above results in the case of Constant Relative Risk Aversion (CRRA) 
utility functions. 
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Example 4.1 The case of CRRA Utility function 

We consider utility functions in the form 



U(x) 



x H 



P 



p<l,p^0, x > 0. 



In this case, we easily compute the optimal wealth process in (|4.2|) : 

-q 



x 



E 



a T {9) 



Z T {0) 
a T {0) 



E 



Z T (8) \~ q 



Zt{0) 



0<t<T, 



where q = The optimal value process is then given for all x > by: 



V e d {x) = -. E a T (9) 
p \ l 



Z T {0) 
a T {9) 



?6 



1-p 



, 6e[0,T]. 



(4.6) 



Notice that the case of logarithmic utility function: U(x) = lnx, x > 0, can be either 
computed directly, or derived as the limiting case of power utility function case: U (x) = 



x 1 as p goes to zero. The optimal wealth process is given by: 



x 



E[a T (e)\FoY Zt(0) 
and the optimal value process for all x > 0, is equal to: 

V e d (x) 



<t<T, 



E[a T 



(0)1^] In ( 



+ E 



6 G [0,T]. 



4.2 The global before-default optimization problem 



In this paragraph, we focus on the resolution of the optimization problem (|3.9|) . We already 
know the existence of an optimal strategy 7r F to this problem, see Remark 13. 1\ and our main 
concern is to provide an explicit characterization of the optimal control. 

We use a dynamic programming approach. For any t £ [0,T], u G Aw, let us consider 
the set of controls coinciding with v until time t: 



K G Aw : 7T F M = u At }. 



Under the standing condition that Vo < co, we then introduce the dynamic version of the 
optimization problem (|3.9p by considering the family of F-adapted processes: 



Vt(u) 



ess sup E 



U{X ¥ T )G T + v e d {x ¥ e (i-4 le ))de 



0<t<T, 



so that Vo = Vq{u) for any v G Aw- In the above expression, X ¥ is the wealth process of 
dynamics (13. 4j) . controlled by 7r F G A(t,v), and starting from Xq. We also denote X V ' V 
the wealth process of dynamics (|3.4|) . controlled by G ^4f 5 starting from Xo, so that it 
coincides with X v until time t, i.e. = X\ t . From the dynamic programming principle 
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(see El Karoui [2]), the process {Vt(v), < t <T} can be chosen in its cad-lag version, and 
is such that for any v £ Aw- 

{Vt{v)+ [ Vf(X£ w (l - vele))d9, < t < T) is a (P, F) — supermartingale. (4.7) 
Jo 

Moreover, the optimal strategy 7r F to problem Vo, is characterized by the martingale pro- 
perty: 

{V t (7r ¥ ) + / V 9 d (xf ' F (l - 7Tele))d6, 0<t<T} is a (P,F) - martingale. (4.8) 
Jo 

In the sequel, we shall exploit this dynamic programming properties in the particular 
important case of constant relative risk aversion (CRRA) utility functions. We then consider 
utility functions in the form 

X P 

U(x) = — , p < 1, p ^ 0, x > 0, 
p 

and we set q = j^-. Notice that we deal with the relevant economic case when p < 0, i.e. 
the degree of risk aversion 1 — p is strictly larger than 1. This will induce some additional 
technical difficulties with respect to the case p > 0. For CRRA utility function, V d {x) is 
also of the same power type, see (|4.6p : 



V e d (x) = U{x)K p e with Kg = E a T (6) 



Zt{0) \ 
a T (6)J 



and we assume that Kg is finite a.s. for all 9 G [0, T\. The value of the optimization problem 
()3.9p is written as 

V Q = sup E[U(xf)G T + [ U(X v /)(l-v e ie) p K p g de], 
ueA v Jo 

In the above equality, we may without loss of generality take supremum over Aw(U), the 
set of elements v € Aw such that: 

E[U(xf)G T + ! U{Xf)(l-u ele yKlde] > -co, (4.9) 
Jo 

and by misuse of notation, we write Av = Aw(U). For any v E Aw with corresponding 
wealth process X u ' governed by (|3.4|) with control v, and starting from Xq, we notice that 
the cad-lag F-adapted process defined by: 



Vtiy) 

[/(At 



* ■= T^k (4-10) 



p ess sup E 



A f Jt A. t 



0<t<T 



does not depend on v € ^4f- It lies in the set L + (¥) of nonnegative cad-lag F-adapted 
processes. Let us also denote by Lf oc (W) the set of F-adapted process 4> s.t. J \4> t \ 2 dt < 
oo a.s. 

We have the following preliminary properties on this process Y. 
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Lemma 4.1 The process Y in (|4.10p is strictly positive, i.e. P[Yj > 0, < t < T] 
Moreover, for all v G Aw, the process 



1. 



&0O 



U{xf)Y t +fu{xf){l-vele) p K p e d6, 0<t<T, (4.11) 
Jo 



is bounded from below by a martingale. 

Proof. (1) We first consider the case p > 0. Then, 



Y, 



ess sup E 



X 



u.¥ 



G T + 



T 



xl 



X 



v,¥ 



(i-4 le fK p e de 



(4.12) 



> E[G T + J K p e d6 T t ] > 0, V*g[0,T], 



by taking in (|4.12j) the control process 7r F G Aw(t, u) defined by 7r F = u s l s <f Moreover, 
since U(x) is nonnegative, the process £ V (Y) is nonnegative, hence trivially bounded from 
below by the zero martingale. 
(2) We next consider the case p < 0. Then, 



Y 



ess inf E 



y¥ 



> Jf.-- 



ess inf E 



G T + 

y¥ 



xl 



X 



x 



Ft 



{\-4leYK p 9 dB 
Vt G [0,T]. 



(4.13) 



Notice that the process J can be chosen in its cad-lag modification. Let us show that for 
any t G [0, T], the infimum in Jt is attained. Fix t G [0, T], and consider, by a measurable 
selection argument, a minimizing sequence (vr n ) n G Aw(t, v) to Jt, i.e. 



lim E 



XX 



X 



v.¥ 



Gt 



Jt 



a.s. 



(4.14) 



Here X n denotes the wealth process of dynamics (|3.4p governed by n n . Consider the (local) 
martingale density process 



exp 



u F 1 
at U 2 



U ¥ 2 



0~Z 



du 



t<s<T. 



By definition of Z l and Ito's formula, the process {Zj,X™,t < s < T} is a nonnegative 
fm , (J 7 s )t<s<T)-loc&l martingale, hence a supermartingale, and so ^[X^Z^Tt] < X^Z\ = 



X^' F . By Komlos Lemma applied to the sequence of nonnegative J-r-measurable random 



variable (X^) n , there exists a convex combination Xj\ G conv(^^ 1 , k ^ tx} such that (X?p^ n 
converges a.s. to some nonnegative J-r-measurable random variable Xt- By Fatou's lemma, 
we have Xt '■= EfX^Zyl^] < -X^' ■ Moreover, by convexity of x — > x p , and Fatou's lemma, 
it follows from (|4.14p that 



Jt > E 



Xj 1 

xF 



Gt 



Ft 



a.s. 



(4.15) 



Now, since p < 0, Jt < oo, and Gt > a.s., we deduce that Xt > 0, and so X t > a.s. 



Consider the process X* 



x t 



■E[=£X T \F S ], t < s <T. Then, < s < T} is a 
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strictly positive (P, (,F s )t<< i <T)-niartingale, and by the martingale representation theorem 
for brownian filtration, using same arguments as in the end of proof of Theorem 14.11 we 
obtain the existence of an (J r s ) t < s <T-&d&pted process 7f* = (7f*)t< s <T satisfying J* f T |w*crf \ 2 ds 
< oo, such that X* satisfies the wealth process dynamics (j3.4|) with portfolio 7f on (t,T), 
and starting from X t = X^' F . By considering the portfolio process it G „4f(^^) defined 



by ir s = u s l s <.t + ^l^s>t, for < s < T, and denoting by X n,¥ the corresponding wealth 



X^i 



process, it follows that Xj' = X* for t < s < T, and in particular Xj' 
Xy a.s. From (|4.15|) . the nonincreasing property of x — > x p , and definition of Jf, we deduce 



that 



J/ 



E 



X, 



7f,F. 

T 



a.s. 



(4.16) 



and as a byproduct that X^' F = X7 1 . The equality (|4,16p means that the process J = 
(Jt)te[o,T] 1S a modification of the process J = (Jt)te[0,T]- Since, J and J are cad-lag, they 
are then indistinguishable, i.e. P[Jt = Jj, < t < T] = 1. We deduce that the process J, 
and consequently Y, inherit the strict positivity of the process J. 
From (l4~T3l) . we have for all v E Aw, t e [0, T], 



ercn 



ess sup E 



U{Xl)G T + 



fu(xl)( 

Jo 



l-4 le ) p K p e dd 



> M?:=E U{X V /)G T + I U{X v /){l-u ele ) p K p e de 



Ft 



(4.17) 
t G [0,71, 



by taking in (|4.17p the control process ir ¥ = v G A&(t, v). The negative process (M t ") ig r 0j r] 
is an integrable (recall (|4.9p ) martingale, and the assertions of the Lemma are proved. □ 

In the sequel, we denote by L b + {¥) the set of processes Y in L + (F), such that for all 
v G Ay, the process £ U (Y) is bounded from below by a martingale. The next result gives 
a characterization of the process Y in terms of backward stochastic differential equation 
(BSDE) and of the optimal strategy to problem (|3.9p . 



Theorem 4.2 When p > (resp. p < 0), the process Y in (|4.10p is the smallest (resp. 
largest) solution in L b + {¥) to the BSDE: 



T 



G T + / f(e,Y e , 



Y, 



for some <j) G L 2 oc (W), and where 



f(t,Y t ,, 



p ess sup 



fadWe, 0<t<T, 



VY t W t a!\ 2 + K? { ±- 



vtlt) 



(4.18) 



(4.19) 



The optimal strategy ("ftt)te[o,T] t° problem ()3.9p attains the supremum in (|4. 19|) . Moreover, 



under the integrability condition: f Q 
taken pointwise, i.e. 



T 



Ki 



2p_ 
2-p 



dt < oo a.s., the supremum in (|4.19p can be 



p ess sup 

n<l/"/t 



^Y^f + Kf^ 



Kit) 
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while the optimal strategy is given by: 



71 + 



arg max 

T<l/7t 



{fj, t Y t + a t (j)t)TT — l^|vrcr t | + iff 



and satisfies the estimates: 



it; 



(i- P )y t |of|2 



where 



7T, 



mm 



i 

2-p 



(i-p)kfl 2 



< 7T 4 < 7T 4 , 



+ 



P 



<t <T, 



0<t<T, 



(4.20) 



(i-p)*W 



7i 



Proof. By Lemma |4. H the process V lies in L^(F). From (|4~T|) . we know that for any v G 
Aw, the process ^(y) is a (P, F)-supermartingale. In particular, by taking v = 0, we see that 
the process {y + J Kgd9,0 < t < T} is a (P, F)-supermartingale. From the Doob-Meyer 
decomposition, and the (local) martingale representation theorem for brownian motion 
filtration, we get the existence of <f> G L? C (W), and a finite variation F-adapted process A 
such that: 



dY t 



dW t - dA t , < t < T. 



(4.21) 



From (|3.4p and Ito's formula, we deduce that the finite variation process in the decompo- 
sition of the (P, F)-supermartingale £"(y), v G At, is given by —A u with 



dAl 



-dA t 



P 



1 — p 



Y t \v t a t \ 2 + K\ 



p 



dt 



Now, by the supermartingale property of £ U (Y), v G Aw, which means that A v is nonde- 
creasing, and the martingale property of £ n¥ (Y), i.e. A = 0, this implies: 



dA t 



P 



^Y + af^nj 



1 — p 



V l^ F ^ F l 2 _i_ jsP 



(1 - *ht) p 



dt 



p ess sup 

ueAw 



(^Y + aUt^t 



1 — p 



Y t \v t o?\ 2 + K 



(i - u t j t y 



p 



dt. 



Observing from (|4.10p that Yr = Gt, this proves together with (|4.21|) that (Y, (p) solves the 
BSDE (|4.18p . In particular, the process Y is continuous. 

Consider now another solution (Y, <f>) G L b + (¥) x L 2 loc (W) to the BSDE ([4TTg]) . and define 
the family of nonnegative F-adapted processes (, U (Y), v G Av, by: 



U{xf)Y t + I U{X v /){l-v ele fKlde, 0<t<T. 



(4.22) 



By Ito's formula, we see by the same calculations as above that: d^{Y) = dM" — dA v t , 
where A v is a nondecreasing F-adapted process, and M v is a local (P, F)-martingale as a 
stochastic integral with respect to the brownian motion W . By Fatou's lemma under the 
condition Y G L5_(F), this implies that the process £, U (Y) is a (P, F)-supermartingale, for 
any v G Aw- Recalling that Yr = Gt, we deduce that for all v G Aw 



E 



U(xf)G T + I U{X u /){l-veleYKde F t < U{X v /)% < i < T. 



T 
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K t 


f 


of 



If p > (resp. p < 0), then by dividing the above inequalities by U (X^ ,¥ ), which is positive 
(resp. negative), we deduce by definition of Y (see (|4.12|) and (|4.13|) ). and arbitrariness 
of v € Af, that Yt < (resp. >) 1^, < t < T. This shows that Y is the smallest (resp. 
largest) solution to the BSDE (gjgj) . 

Next, we make the additional integrability condition: 

2 ~ p cft < oo, o.s. (4.23) 

Let us consider the function F defined on {(uj,t, it) € x [0, T] x R : 7r < 1/7^(0;)} by: 

z P 
(As usual, we omit the dependence of i 7 inw). By definition, we clearly have almost surely 

-f(t,Y t ,<pt) < esssupF(t,vr), < t < T. (4.24) 

P 7T<l/7i 

Let us prove the converse inequality. Observe that, almost surely, for all t 6 [0,T], the 
function it — > F(t,ir) is strictly concave (recall that the process Y is strictly positive), C 2 
on (—00,1/74), with: 

— (t,7T) = ( f 4Y t +af ( f>t)-(l-p)Y t \af\ 2 7T- lt Kf(l-7T 7t r- 1 , 
and satisfies: 

dF OF 
lim F{t,ir) = —00, lim — — (t,vr) = 00, lim — — (t, 7r) = —00. 

tt-^-oo 7T^-00 OTT 7T— >l/7t 07T 

We deduce that almost surely, for all t £ [0, T], the function 7r — > F(t,ir) attains its maxi- 
mum at some point ttJ, which satisfies: 

gfttf - 0. 

By a measurable selection argument, this defines an F-adapted process tt ¥ = (ttJ )te[o,T]- I n 
order to prove the equality in (|4.24p . it suffices to show that such 7r F lies in ^4f, and this 
will be checked under the condition (|4.23p . For this, consider the F-adapted processes tt m 
and tt m defined by: 



7T 



A*t fit M ■ (~M 1 N 

(l-p)|crf|^ (l-p)Y t af Jt 



When TT^ioo) < 1/72(0/), we have: 

and so by strict concavity of F(t, n) in it: ttJ < tt^ . When tt^ 1 (u) > 1/7^(0;), and since ttJ 
< l/7t, we get: jtj < fr^ . Consequently, we have the upperbound: 7rf < ttJ^ , for all £ € 
[0, T]. Notice that by (|2.5p . continuity of the path of Y, and since (p € Lf oc (W), we have: 
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Jo \^t^ a t \ 2 dt < oo a.s. Moreover, since jt > 0, we have |-7r A/ | < |vr A/ |, and thus ir M lies in 
Aw- Next, consider the F- adapted process tt defined by: 

7f t = n^-pt, 0<t<T, 

for some F-adapted nonnegative process p = (pt)te[o,T] to be determined. When jrf 1 (oj) < 
1/74(0;), we have 

r) F 

^(i,7f t ) = (l- p )Y t \af\ 2 p t - lt Kf(l-^ I lt + p tl t) p - 1 

> (l-p)Y t \aJ\ 2 Pt - lt Kf(pat) p - 1 . (4.25) 

When tt^(lu) > l/7t(u;), the inequality (|4.25p also holds true. Hence, by choosing p such 
that the r.h.s. of (|4,25p vanishes, i.e. 

[l-p)Y t \af 



Pt = ( T-, F| 2 ) > o<t<r, 



dF 

we obtain almost surely: -7— (i, 7Tj) > 0, < t < T, and so by strict concavity of F in tt: ixt 

OTT 

< frf. Finally, under (|4.23p . and recalling that Y is continuous, 7 < 1, we easily see that p 
satisfies the integrability condition: \ptaf\ 2 dt < 00 a.s., and so tt lies in Aw- Therefore, 
we have proved that tt ¥ lies in Aw, and satisfies the estimates (|4.20p . □ 

Remark 4.2 The driver f(t,Yt,4>t) of the BSDE (I4.18P is in general not Lipschitz in the 
arguments in (Y~t, (pt), and we are not able to prove by standard arguments that there exists 
a unique solution to this BSDE. 

Remark 4.3 We make some comments and interpretation on the form of the optimal 
before-default strategy. Let us consider a default-free stock market model with drift and 
volatility coefficients /i F and <r F , and an investor with CRRA utility function U(x) = x p /p, 
looking for the optimal investment problem: 

V M = sup E[U(X%)], 

where X v is the wealth process in (I3.4p . In this context, Aw, defined in (13. 3p is interpreted 
as the set of trading strategies that are constrained to be upper-bounded (in proportion) 
by 1/7*. In other words, Vq 1 is the Merton optimal investment problem under constrained 
strategies. By considering, similarly as in (|4.10p . the process 



— V 



Y t = p ess sup E 



U 



X 



xf 



0<t<T, 



and arguing similarly as in Theorem 14.21 one can prove that Y M is the smallest solution to 
the BSDE: 



Y t M = 1 + J f M (0,Y e M ,4 J )de- J t 4'dWe, 0<t<T, 
for some 4> M £ Lf oc (W), where 

f M {t,Yf,tf) = pesssupUpjY^ + af^-^Y^lnafl 2 
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while the optimal strategy for V^ 1 is given by: 



7T* = mm — ^ + 



Notice that when the coefficients fi ¥ , a ¥ and 7 are deterministic, then Y M is also deter- 
ministic, i.e. <p M = 0, and is the positive solution to the ordinary differential equation: 



Y t M = 1 + J f M (6,Y 9 M )d6, 0<t<T, 



with f M (t,y) = pysup JT<1/lt [^7r - ^{naf | 2 ] =: pyc(t), and so: Y t M = exp(p c(9)d6). 

F 

Moreover, the optimal strategy is irf 1 = min ( p^yHrp > m particular, when there 

is no constraint on trading strategies, i.e. 7 = 0, we recover the usual expression of the 



,7 



optimal Merton trading strategy: tt^' 1 = ^^^rp ■ 

Here, in our default stock market model, the optimal before-default strategy tt v satisfies 
the estimates (I4.20j) . which have the following interpretation. The process ir M has a similar 
form as the optimal Merton trading strategy tt m described above, but includes further 
through the process Y and K, the eventuality of a default of the stock price, inducing a 
drop of size 7, and then a switch of the coefficients of the stock price from (// F ,<r F ) to 
(fi d ,a d ). The optimal trading strategy 7r F is upper-bounded by ir M , and when the jump 
size 7 goes to zero, it converges to tt m , as expected since in this case the model behaves as 
a no-default market. 

4.3 Example and numerical illustrations 

We consider a special case where fi ¥ , er F , 7 are constants, [i d {0) (J d {9) are only deterministic 
functions of 9, and the default time r is independent of F, so that ctt{9) is simply a known 
deterministic function a{9) of 9 £ R+, and the survival probability G(t) = P[r > t\Tt] = 
P[t > t] = J t °° a(9)d9 is a deterministic function. We also choose a CRRA utility function 
U(x) = p < 1, p ^ 0, x > 0. Notice that V^s) = t/^M) = U(x)k(6) p with 



a r (0) 



= (Ea T (0)(^ZT) " ) ? =a(0)^exp(^ -— (T 







2 


(i 




1 



V 



Moreover, the optimal wealth process after-default does not depend on the default time 
density, and the optimal strategy after-default is given, similarly as in the Merton case, by: 



(i-p)k"(»)l 



2 ' 



On the other hand, from the above results and discussion, we know that in this Markovian 
case, the value function of the global before-default optimization problem is in the form Vo 
= v(0,X ) with: 

v{t,x) = U(x)Y(t), 

where Y is a deterministic function of time, solution to the first-order ordinary differential 
equation (ODE): 

Y(t) = G(T) + J T f(0,Y(O))dB, te[0,T], (4.26) 
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with 



7T<l/7 



f(t,y) = psup l/K-^-P\na v \ 2 )y + k(tT { - (4.27) 



P 



There is no explicit solutions to this ODE, and we shall give some numerical illustrations. 

The following numerical results are based on the model parameters described below. 
We suppose that the survival probability follows the exponential distribution with constant 
default intensity, i.e. G(t) = e~ xt where A > 0, and thus the density function is a(8) = 
Ae~ Ae . The functions [i d {9) and a d (9) are supposed to be in the form 

Ao) = Ao) = ^ F (2-|), ee[o,T), 

which have the following economic interpretation. The ratio between the after and before- 
default rate of return is smaller than one, and increases linearly with the default time: the 
after-default rate of return drops to zero, when the default time occurs near the initial 
date, and converges to the before-default rate of return, when the default time occurs near 
the finite investment horizon. We have a similar interpreation for the volatility but with 
symmetric relation: the ratio between the after and before-default volatility is larger than 
one, decreases linearly with the default time, converging to the double (resp. initial) value 
of the before-default volatility, when the default time goes to the initial (resp. terminal 
horizon) time. 

To solve numerically the ODE (|4.26p . we apply the Howard algorithm, which consists 
in iterating in (j4.27j) the control value tt at each step of the ODE resolution. We initialize 
the algorithm by choosing the constrained Merton strategy 

£lM .,( n" i 



7r = mm 



(l-p)|a F | 2 ' 7 



In the following Table [U we show the impact of the loss given default 7 on the optimal 
strategy ttJ . The numerical tests show that except in some extreme cases where both 
the default probability and the loss given default are large, the optimal strategy is quite 
invariant with respect to time t in most cases we consider. So we give below the optimal 
strategy as its expected value on time. We perform numerical results for various degrees of 
risk aversion 1 — p: smaller, close to and larger than one, and with /jF = 0.03, a ¥ = 0.1, T 
= 1 and A = 0.01. 



Table 1: Optimal strategy vs constrained Merton. 





P = 


0.2 


P ~ 


■+ 


P = 


-0.2 


7 








7T M 






0.01 


3.73 


3.74 


2.99 


3.00 


2.49 


2.50 


0.1 


3.57 


3.74 


2.86 


3.00 


2.38 


2.50 


0.5 


1.58 


2.00 


1.38 


2.00 


1.22 


2.00 


0.8 


0.91 


1.25 


0.80 


1.25 


0.70 


1.25 



First, observe that, when we take into account the counterparty risk, the proportion 
invested in the stock is always smaller than the Merton strategy without conterparty risk. 
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Table 2: Optimal strategy with various A and 7. 









P = 


0.2 


P ~ 


-> 


P = 


-0.2 








7 = 0.1 


7 = 0.5 


7 = 0.1 


7 = 0.5 


7 = 0.1 


7 = 0.5 


7T M 


3.75 


2.00 


3.00 


2.00 


2.50 


2.00 


A = 0.01 


PD= 


0.01 


3.57 


1.58 


2.86 


1.38 


2.38 


1.22 


A = 0.05 


PD= 


0.05 


2.93 


0.26 


2.35 


0.21 


1.96 


0.18 


A = 0.1 


PD= 


0.10 


2.22 


-0.90 


1.78 


-0.70 


1.49 


-0.58 


A = 0.3 


PD= 


0.26 


0.00 


-3.96 


0.00 


-3.00 


0.00 


-2.40 



Secondly, the strategy is decreasing with respect to 7, which means that one should reduce 
the stock investment when the loss given default increases. 

Next, we examine the role played by the default intensity A. In Table [2j the column PD 
represents the default probability of the counterparty up to T with the given intensity A, i.e. 
PD = P(r < T) = 1 — e~ XT . We see that the stock investment decreases rapidly when the 
default probability increases. Moreover, when both default probability of the counterparty 
and the loss given default of the stock are large, one should take short position on the stock 
in the portfolio investment strategy before the default of the counterparty. Then at the 
default time r = 6, the optimal strategy is switched to 7t d (6), which is always positive. 

We also compare the value function obtained in our example to that in the classical 
Merton model, that is, the solution Y(t) to the ODE (|4.26p and the function Y M (t) deduced 
with k(t) = and G(T) = 1. In Figure [H the curves represent different values of 7 such 
that tt m = nz^ynrp < y- The value function Y(t) obtained with counterparty risk is 
always below the Merton one. It is decreasing on time and also decreasing w.r.t. the 
proportional loss 7. In addition, for a given default intensity A, all curves converge at T to 
G(T) = e- XT . In Figure El the curves represent different values of A. The value function 
Y(t) is also decreasing w.r.t. the default intensity A. However, the final value of each curve 
corresponds to different values of Gt- 
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Figure 1: Value function - optimal vs Merton: p = 0.1, A = 0.01. 
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Figure 2: Value function - optimal vs Merton: p = 0.1, 7 = 0.1. 
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5 Conclusion 



This paper studies an optimal investment problem under the presence of counterparty risk 
for the trading stock. By adopting a conditional density approach for the default time, 
we derive a suitable decomposition of the initial utility maximization problem into an 
after-default and a global default one, the solution to the latter depending on the former. 
This makes the resolution of the optimization problem more explicit, and provides a fine 
understanding of the optimal trading strategy emphasizing the impact of default time 
and loss given default. The density approach can be used for studying other optimal 
portfolio problems, like the pricing by indifference-utility, with counterparty risk. A further 
topic is the optimal investment problem with two assets (names) exposed both to bilateral 
counterparty risk, and the conditional density approach should be relevant for such study 
planned for future research. 
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